Velocity dependent interactions and a new sum rule in bcc He 
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Recent neutron scattering experiments [PRL, 88, p. 195301 (2002)] on solid 4 He, discovered a new 
optic-like mode in the bcc phase. This excitation was predicted by a recently proposed model that 
describes the correlated atomic zero-point motion in bcc Helium in terms of dynamic electric dipole 
moments. Modulations of the relative phase of these dipoles between different atoms describes the 
anomalously soft Ti(110) phonon and two new optic-like modes, one of which was recently found in 
the neutron scattering experiments. In this work we show that the correlated dipolar interactions 
can be written as a velocity dependent interaction. This then results in a modified f-sum rule for 
the Ti(110) phonon, in good agreement with the recent experimental data. 
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The phonons in a quantum solid, such as solid He, show 
various quantum effects arising from the unusually large 
zero-point motion of the atoms. The strong quantum 
pressure in solid He lowers the density and introduces 
large anharmonic terms in the interatomic potential. In 
bcc He the potential even has a double- well feature |l|, |) . 
Describing the phonon spectra of this highly anharmonic 
quantum solid in a consistent manner has therefore pre- 
sented a theoretical challenge over the years. 

The quantum effects appear in the first moment of the 
dynamic structure factor S(to,k), as measured by neu- 
tron scattering ||. The Debye- Waller factor d(k) is un- 
usually small due to the large spread in the atomic po- 
sitions. In addition there have been measurements || 
of apparent oscillations of the Debye- Waller factor as a 
function of the scattering wavevector k. These oscilla- 
tions seemed to be confined to some modes, while ab- 
sent in the transverse Ti(110) mode. The explanation 
of these unusual oscillations was given in terms of inter- 
ference effects between the single and multi-phonon ex- 
citations, whereby intensity is transferred between them 
[jlj . The mixing of the phonon modes is attributed to the 
anharmonicity of the interatomic potential. From sym- 
metry considerations these types of calculations predict 
that the oscillatory contribution to d{k) vanishes both at 
the Brillouin zone center and boundary (z.b.). 

Recently there has been a new neutron scattering ex- 
periment Q in solid 4 He. The main surprise in this ex- 
periment was the discovery of a new, gapped excitation 
branch in bcc 4 He. This excitation was found at the en- 
ergy predicted by a new model of the dynamic effects 
of the zero-point motion on the phonon spectra ^). 
Usually a bcc lattice is not expected to support optic- 
like phonons, so the measurement of this new excitation 
branch provides support for the basic idea of the new 
model, namely the introduction of an additional degree 
of freedom in bcc He. Within the new model, this de- 
gree of freedom is the relative phase of a correlated part 
of the zero-point motion between the atoms. The latest 
experiment has also provided data which shows appar- 
ent oscillations of the Debye- Waller for the transverse 



Ti(110) mode §. It is this data which we shall analyze 
in this letter, in terms of the new model. 

In previous papers ^ ^ we have presented a new 
model to describe some of the correlations in the atomic 
zero-point motion in bcc Helium. The main point is that 
there is an anharmonic, low energy mode of atomic mo- 
tion along the major axes of the bcc lattice, which is not 
described by the harmonic approximation. In bcc He 
this mode has energy Eq ~ 3 — 10K Q, and is the only 
free parameter in our model, which we take from empir- 
ical data. Additionally, the low mass of the He atom al- 
lows relatively large breaking of the Born-Oppcnheimer 
approximation which results in motion-induced electric 
dipoles, due to relative displacement of the nuclei with 
respect to the electronic cloud ^ ^ . When these dipoles 
are uncorrelated they introduce a negligible addition to 
the usual Van-der Waals interaction. By contrast, cor- 
related and highly directional atomic motion results in 
dipolar interactions between these dipoles, which is of 
the order of the energy of the local atomic motion Eq. 
There can therefore arise a state of quantum resonance 
between these interacting electric dipoles and the atomic 
motion. A resonant state between the atomic motion and 
electric dipolar interactions lowers the ground-state en- 
ergy of the system. In the ground-state of the bcc phase 
we found g that the dipolar interaction is minimized 
when the dipoles have the phase relation shown in Fig.l. 

The dipolar interactions between the motion- 
correlated dipoles is pf 
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where the electric dipole moment in bcc He turns out to 
be ||: fj, ~ e • 0.01, and r i0 = r — and the index 
i runs over all atoms. The resonance condition means 
that Eq is the frequency of both the oscillating dipoles 
and the atomic motion, i.e. Edip = ~Eq/2 (~ — 3K in 
bcc 4 He). Note that the dipole moments are dynamic, 
and the expectation value of the dipole moment of the 
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crystal and of each atom is zero, as required by time 
reversal symmetry. It is the matrix element of correlated 
dipoles (poPi) which is non-zero, and appears in ([!]). 

The modulations of the relative phases of the dipoles 
on different atoms coincide with only one phonon mode, 
due to the lower symmetry of the dipolar ground-state 
compared with the bcc lattice (Fig. 1) Q. This mode 
turns out to be the anomalously soft Ti(llO) phonon, 
which agrees very well with the experimental data [Q. 
In addition our model predicts two new modes which 
are gapped (optic- like) [|[ [|. One of these modes was 
recently observed in neutron scattering experiments [Q. 
Since the dipolar degree of freedom involves atomic mo- 
tion, it couples to the atomic density operator. A flipping 
of a dipole with respect to the ground-state configuration 
(Fig.l) results in a dynamic density fluctuation. A flip- 
ping of two adjacent dipoles forms the localized excitation 
of energy 2Eq |2j recently found by neutron scattering [Q . 
Our model therefore gives a good description of some ef- 
fects of the anharmonic potential in bcc He, as they are 
manifested in the excitation spectrum. 

Since the dipolar interactions arise from correlated 
atomic motion, they can be naturally written as a 
velocity-dependent potential. Such a term will mod- 
ify the Nozieres-Pines f-sum rule, which results from a 
Hamiltonian where the momentum operator appears only 
in the kinetic energy term H (i.e. no velocity-dependent 
interactions). In this paper we calculate the f-sum rule 
for the Ti(llO) phonon, which we describe as a phase 
modulation of the correlated dipolar array |||. 

Before proceeding let us make a short comment con- 
cerning the appearance of velocity-dependent interac- 
tions in physics. One example is superfluid 4 He || . In the 
superfluid a velocity-dependent interaction is introduced 
into the energy density functional in order to take into ac- 
count short-range atomic motion correlations. These are 
termed "backflow" and introduce a non-local velocity- 
dependent term to the potential energy. The result is a 
modified f-sum rule for the phonon-roton dispersion in 
bulk superfluid 4 He (Eq.(44) in ||]). All this is very sim- 
ilar to our results, only that we do not introduce any 
unknown fitting parameters, as is done for the superfluid 
case. 

The dipolar interaction matrix element A"(k) in the 
presence of a phase modulation (i.e. a Ti(llO) phonon), 
is the Fourier transform of ([!]) Q 
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shown in Fig. 2. 

By analogy with an oscillating electromagnetic field 
interacting with an atom |To|| , we can write the electric 
dipole moment in term of a local momentum operator 



p = eX -> P ■ 
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where ujq = EqJH is the resonance frequency of atomic 
and dipolar vibration, and the ""symbol is used for oper- 
ators. 

The space representation of the dipolar interaction en- 
ergy (Q) between the zero-point dipoles in the ground- 
state can now be written in terms of the local momentum 
operators 
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where ry = rj — r,j. Note that similar to the treatment 
Q of superfluid 4 He, we find a non-local kinetic term. 

The contribution of an interaction X (Q) to the f-sum 
rule is found by performing the usual double commutator 
with the density operator B 
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where = e %k ' Ti is the standard density operator 
and P is a unit vector of the local momentum so that 
we keep track of the phase relation between the corre- 
lated motion of different atoms. The commutator (||) 
is performed at equal times |]ll). When calculating the 
summation in (|^) we have to specify the phase relation 
between the dipoles on different sites. If there is no corre- 
lation between the dipoles, the interaction X{k) is zero, 
and the correlation function of the atomic displacements 
described by these dipoles, i.e. S(ui,k), is equally zero. 
In the bcc phase we have the dipoles correlated in the 
ground-state as shown in Fig.l, with resulting non-zero 
dipolar interactions (Fig. 2). The momentum operators 
in (^J) replace the usual kinetic energy term for the part 
of the atomic motion which is described by the dipolar 
degree of freedom. 

Using this double commutator (|^) in the usual defini- 
tion of the dynamic structure- fact or S(u>,k) || we 
get that the new sum rule for the Ti(llO) phonon is 
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where just as for usual phonons, the single-phonon contri- 
bution to the f-sum rule is modified by the Debye- Waller 

factor d{k) = exy-^)^ (ACB sum rule)[0, [if], with 
(u 2 ) being the isotropic zero-point spread of the atomic 
wavefunction. We assumed here that the dynamic struc- 
ture factor S(u>,k), i.e. the correlations in atomic dis- 
placements, for the Ti(llO) phonon are completely de- 
scribed by the dipolar degree of freedom alone. This is 
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justified by the excellent agreement between our calcula- 
tion of the and the T^llO) spectrum experimental data 

The new sum rule we find for the Ti (110) phonon (|6j) is 
essentially the usual f-sum rule multiplied by the function 
2\X(k)\/ Eq. We find that at k — > 0, we recover the usual 
f-sum rule, since \X(k — > 0)| — > E /2. While the usual f- 
sum rule follows from particle conservation, the number 
of excited dipoles out of the correlated ground-state is 
not strictly conserved. In particular, when the function 
X(k) goes to zero at the Brillouin z.b., the commutator 
(|J) vanishes. The vanishing of M\{k) can be understood 
as a vanishing of both the restoring force of the dipolar 
interactions and therefore the dipole-induced correlations 
between the atomic motion on different sites. Exciting a 
dipole does not affect the others at this particular mo- 
mentum, and the neutron can not transfer any energy to 
this degree of freedom. Note that M\ < Nh 2 k 2 /2m since 
2|X(A;)| < E , so that the energy constraint which speci- 
fies that the maximum energy transferred by the neutron 
to the atom is the given by recoil energy, is satisfied. 

Before comparing to the experimental data we note 
that in the experiments the neutron scattering is usu- 
ally done in the second or higher Brillouin zone, so the 
periodic behavior of X(k) (Fig. 2) has to be taken into 
account. The experimental set-up is described elsewhere 
H 0], and both sets of data are for bec 4 He. Our calcula- 
tion is equally applicable for bec 3 He 0, for which there 
is at present no neutron scattering data. 

The calculated M\{k) (||) is shown in Fig. 3 compared 
to the neutron scattering data || |). For d(k) we use 
the value (it 2 ) = 0.17A 2 , taken from the measured 
Q Debye temperature of ~ 27K. By plotting the ra- 
tio Mi(k)/ ([k • e(k)] 2 d(k) 2 ) we can easily see deviations 
from the usual f-sum rule, which gives a constant value 
(e(fc) being the phonon polarization) jl). The data in 
Fig. 3 is plotted against the full scattering momentum 
of the Ti(100) phonon as measured in the most recent 
experiments j|, |), so both the older data || and the 
calculation are mapped accordingly. 

We find that our calculation is in good agreement with 
the experimental data. The calculated oscillation in the 
intensity is observed, in contrast with the usual f-sum 
rule result. The observed periodicity is different from the 
one calculated previously for interference terms jlj , which 
gives a correction that vanishes both at the z.b. and the 
zone center. Additionally, the multi-phonon interference 



calculation does not predict the complete vanishing of 
the intensity at any momentum, unlike the result of the 
model presented here. Note that both our model and pre- 
vious calculations attribute the oscillations in the scat- 
tered intensity to anharmonic effects. 

The prediction that the intensity should vanish at the 
Brillouin z.b. is not observed (Fig. 3), only a marked re- 
duction in the scattering intensity. This can arise from 
the finite resolution of the experiment B allowing in 
stray scattering from other phonons with similar energy 
(at the z.b. the Ti(110) branch meets other phonon 
branches). At the z.b. with the smaller momentum 
(k ~ 1.85A -1 ) there is higher scattering intensity then at 
the higher momentum z.b. (k — > 2.85A -1 ), as expected 
for stray scattering from other phonons. Future high res- 
olution data at the zone boundaries is needed to resolve 
this question. 

In this work we presented a modified f-sum rule for the 
anomalously soft Ti(110) phonon mode of bec He. We 
find an apparently oscillating Debye- Waller factor and a 
complete vanishing of the first moment at the Brillioun 
zone boundary. The new f-sum rule is qualitatively dif- 
ferent from the standard response of a harmonic phonon. 
This unusual result is due to correlations in the atomic 
zero-point motion, which can be written in terms of a 
velocity-dependent interaction. The calculated sum rule 
is then compared with recent neutron scattering data, 
which does show some of the predicted features. Com- 
bined with the recent observation of a new excitation 
mode |Q , the sum rule data presented here provides extra 
support for the proposed model of correlated zero-point 
electric dipoles in the ground-state of bec He ||, ^ . 
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FIG. 1: The two degenerate 'antiferroelectric' dipole arrange- 
ments in the ground-state of the bcc phase. The system is 
in quantum resonance between these two configurations with 
resonance frequency u)q = Eq/H. The arrows show the instan- 
taneous direction of the dipoles along one of the major axes. 
Similar dipolar arrays exist along the orthogonal major axes. 
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FIG. 2: The dipolar interaction matrix element X(k) (Eq.g) 
in the bcc phase of 4 He, along the (110) direction. 
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FIG. 3: The ratio Mi(fc)/ ([k ■ e(k)] 2 d(k) 2 ) given by Eq| 
(dashed line), compared with the experimental data for the 
Ti(llO) phonon: crosses Q and stars pj. The usual form 
for a harmonic phonon is given here by the horizontal dotted 
line. 



